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Abstract 

The  objective  of  this  investigation  is  to  evaluate 
the  transient  effects  of  a  reaction  control  jet  on  the 
aerodynamic  performance  of  a  generic  interceptor 
missile.  Three-dimensional  computations  of  the 
highly  turbulent  flow  field  produced  by  a  pulsed,  lat¬ 
eral-jet  control  thruster  and  the  interaction  of  this  jet 
with  the  supersonic  free  stream  and  missile  bound¬ 
ary  layer  were  completed  for  different  altitudes  and 
thruster  conditions.  A  generic  supersonic  missile- 
interceptor  configuration  consisting  of  a  long,  slen¬ 
der  body  (L/D  =  14.1)  containing  fixed  dorsal  and 
tail  fins  was  used  in  this  study.  Parametric  compu¬ 
tational  fluid  dynamic  (CFD)  solutions  were 
obtained  at  two  altitudes  of  64.5  kft  (1 9.7  km)  and 
115  kft  (35.1  km).  Computation  of  the  flow-field 
behaviors  at  each  altitude  were  completed  for  the 
following  assumptions:  (1)  steady-state  conditions, 
lateral  control  jet  turned  off,  (2)  steady-state  condi¬ 
tions,  lateral  control  jet  turned  on,  (3)  transient  tran¬ 
sient  jet  startup  simulation,  and  (4)  transient  jet 
shutdown  simulation.  A  thermally  perfect  gas 
(y  =  1 .4)  was  assumed  for  the  Mach  number  5  free 
stream  and  the  Mach  number  3  lateral  jet.  Vehicle 
forces  and  moments  were  obtained  for  each  solu¬ 
tion  by  integrating  the  surface  pressures  and  vis¬ 
cous  shear  stresses  computed  on  the  missile  sur¬ 
faces.  These  results  are  applied  to  assess  the  influ¬ 
ence  of  the  jet  interaction  (Jl)  effects  on  the  tran¬ 
sient  aerodynamic  performance  of  the  missile.  The 
analysis  indicates  that  strong  transient  influence  is 
predicted  in  the  integrated  normal  force  and  pitch¬ 
ing  moment.  These  effects  may  be  influenced  by 
the  dorsal  fin  interaction  with  the  jet  interaction  (Jl) 
region. 


Introduction 

To  complete  its  mission  successfully,  a  missile 
defense  interceptor  must  be  highly  maneuverable 
when  it  travels  at  supersonic  or  hypersonic  speeds. 
Quick-response  maneuverability,  especially  during 
the  end-game  phase  of  the  interceptor’s  mission,  is 
achieved  by  a  rapid  airframe  response-time  to  the 
attitude  control  system.  Surface-mounted,  fast¬ 
reacting,  lateral-jet  control  thrusters  issuing  at  large 
angles  relative  to  the  interceptor’s  direction  of  flight 
offer  an  effective  supplement  to  conventional  aero¬ 
dynamic  control  surfaces,  producing  the  required 
response  times  and  improving  the  missile’s  agility 
and  maneuverability. 

The  mutual  interference  of  the  jet  thruster 
exhaust  with  the  supersonic  free  stream  leads  to 
thrust  and  moment  amplifications  due  to  high-pres¬ 
sure  regions  that  form  upstream  of  the  jet  on  the 
missile  surface.  This  high-pressure  region  is  cre¬ 
ated  by  the  shock  structure  that  develops  in  the 
supersonic  free  stream  in  front  of  the  jet.  Large 
regions  of  separated  flow  created  by  the  missile 
boundary-layer  interactions  with  the  shock  cause 
the  high-pressure  area  to  increase  in  size.  This 
effect  amplifies  the  response  of  the  divert  jets. 
Effective  exploitation  of  this  effect  can  lead  to 
improved  effectiveness  of  the  missile  interceptor. 
These  solutions  demonstrate  that  CFD  can  be 
applied  to  acquire  a  more  complete  understanding 
of  this  behavior.  Insight  into  this  phenomenon  is 
necessary  because  at  some  missile  orientations 
and  flow  conditions,  the  mutual  interference  of  the 
jet  thruster  flow  field  with  the  free  stream  leads  to 
“deamplification”  (i.e.,  negative  effects).  Therefore, 
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an  understanding  of  the  controlling  factors  that  pro¬ 
duce  thrust  amplification,  as  well  as  those  that  pro¬ 
duce  thrust  deamplification,  is  critical  to  developing 
a  credible  design  basis  for  optimal  missile  perfor¬ 
mance.  The  CFD  results  from  this  investigation, 
complemented  with  ground  simulations  can  be 
applied  to  understand  and  quantitatively  assess 
this  phenomenon.  The  Jl  phenomena  are 
described  in  Refs.  1  and  2. 

The  CFD  solutions  completed  for  this  investiga¬ 
tion  demonstrate  the  potential  to  numerically  simu¬ 
late  complex,  three-dimensional  Jl  flow  fields  and 
assess  the  influence  of  the  transient  jet  interaction 
effects  on  the  aerodynamic  performance  of  the 
missile.  However,  just  as  Jl  phenomena  pose  a 
challenge  for  ground-test  simulations,  it  is  a  formi¬ 
dable  area  for  numerical  simulations  as  well. 
Numerically  simulating  the  transient,  lateral  control 
jet  thruster  flow-field  behavior  as  it  merges  and 
interacts  with  the  supersonic  free-stream  flow  is 
one  of  the  most  complex  propulsion  flow-field  phe¬ 
nomena  to  be  addressed  using  state-of-the-art 
CFD  techniques.  A  numerical  scheme  applied  to 
transient  Jl  flow-field  problems  must  include  the 
capability  to  accurately  resolve  the  complex  shock 
structures  and  boundary-layer  separation  regions 
associated  with  such  flows,  and  it  also  predicts  the 
time-accurate  characteristics  of  these  highly  vorti¬ 
cal  and  possibly  chemically  reacting  flows.  (This 
study  does  not  include  chemical  effects.) 

Objectives 

The  objectives  of  this  effort  are: 

1 .  To  demonstrate  the  feasibility  of  CFD  tech¬ 
niques  to  analyze  Jl  flow  transient  behaviors  and 
verify  the  utility  of  the  solutions  for  interceptor  mis¬ 
sile  design  and  development  concepts. 

2.  To  provide  an  extensive  interrogation  of 
computational  flow-field  results  obtained  at  repre¬ 
sentative  operational  conditions.  Detailed  analy¬ 
ses  of  transient  Jl  flow-field  behaviors  are  not 
plausible  from  the  limited  ground  test  collections 
measuring  forces  and  moments.  The  transient 
effects  are  limited  by  the  response  of  the  instru¬ 
mentation,  and  spatial  flow-field  visualization  of 
complex,  time-accurate  transient  events  has  not 
been  captured  in  ground  simulations. 


3.  To  computationally  visualize  and  quantify  the 
transient  Jl  effects  resulting  from  a  pulsed,  lateral- 
jet  control  thruster  on  the  aerodynamic  perfor¬ 
mance  of  a  generic  supersonic/hypersonic  ballis¬ 
tic  missile  defense  interceptor  at  representative 
operational  conditions. 

Approach 

The  approach  consisted  of  a  computational 
experiment  of  Jl  flow-field  behaviors  involving 
parametric  CFD  investigations  of  the  computed 
flow  fields  and  resulting  forces  and  moments  at  two 
altitude  conditions.  The  effects  of  the  divert  jet 
under  transient  startup  and  shutdown  are  evalu¬ 
ated,  as  well  as  steady-state  Jl  behaviors.  Steady- 
state  solutions  were  also  obtained  when  the  divert 
thruster  was  present  on  the  missile  body,  but  not 
operating.  The  same  missile  and  divert  thruster 
geometries  were  assumed  for  all  simulations.  The 
generic  missile  geometry  shown  in  Fig.  1  was  pro¬ 
vided  by  the  John  Hopkins  University  Applied 
Physics  Laboratory.  A  generic  long,  slender  missile 
body  design  including  forward  and  aft  fins  was 
evaluated.  This  geometry  is  representative  of  inter¬ 
ceptor  designs  and  was  selected  because  divert 
thruster  transient  effects  are  suspected  to  be 
potentially  influenced  by  the  presence  of  fins  on 
the  missile.  The  missile  free-stream  Mach  number 
was  maintained  at  Mach  number  5.0.  The  Mach 
number  at  the  exit  of  the  diverging  divert  thruster’s 
nozzle  was  nominally  3.  A  thermally  perfect  gas 
was  assumed  in  the  CFD  approximation  for  the 
free  stream  and  the  divert  flow.  The  ratio  of  specific 


Fig.  1 .  Surface  definition  showing  zonal 
decomposition. 
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heats  (y)  equal  to  1.4  was  specified  for  the  free 
stream  and  the  divert  jet  flow.  Three-dimensional, 
thermally  perfect  flow-field  solutions  were  com¬ 
pleted  for  the  interceptor/divert  motor  flow-field 
interactions  at  two  altitude  conditions,  64.5  left 
(19.7  km)  and  115  kft  (35.1  km).  These  missile 
flight  conditions  are  representative  of  the  actual 
operation  scenario  of  the  missile  interceptor.  The 
angle  of  attack  of  the  interceptor  was  specified  as 
0.0  deg.  The  divert  jet  was  maintained  at  the  same 
position  on  the  interceptor  missile  with  its  effluents 
situated  at  90  deg  to  the  thrust  vector  of  the  inter¬ 
ceptor  for  all  simulations. 

The  resulting  CFD  flow-field  solutions  are  con¬ 
trasted  to  determine  the  transient  effects  of  divert 
thruster  jet  startup  and  shutdown  behaviors.  The 
flow  field  surrounding  the  interceptor  missile  body 
when  the  divert  thrusters  are  present  on  the  inter¬ 
ceptor  body,  but  not  operating  is  also  computed. 
Even  when  the  divert  jets  are  turned  off,  there  is  an 
influence  on  the  missile  body  forces  due  to  cavity 
flow  created  by  the  presence  on  the  thruster  open¬ 
ing  on  the  body.  This  analysis  was  completed  at 
two  altitudes  to  determine  the  altitude  effect,  if  any, 
on  the  interaction  behavior  and  resulting  missile 
body  forces  and  moments. 

The  divert  jet  is  located  between  two  of  the  dor¬ 
sal  fins  at  the  center  of  gravity  for  the  interceptor 
missile.  The  divert  jet  thrust  chamber  consisted  of 
a  converging/  diverging  nozzle  geometry  and  stag¬ 
nation  conditions  specified  such  that  the  nominal 
Mach  number  at  the  nozzle  exit  is  3.0.  The  divert 
jet  flow  was  included  in  the  computational  domain. 
The  stagnation  conditions  specified  for  the  divert 
thrust  are  given  below: 

Poj  =  2590  psia  (1 76  atm) 

Toj  =  41 00°R  (2278  K) 

The  free-stream  conditions  for  the  two  altitudes 
investigated  at  19.7  km  and  35.1  km  are  shown  in 
Table  1 . 

The  CFD  solutions  were  designed  to  study  the 
flow-field  behaviors  focusing  on  the  transient 
effects  on  the  forces  and  moments  that  result  from 
the  operation  of  a  Mach  3.0  pulsed  jet,  located  at 


Table  1.  Free-stream  Conditions 


Parameter 

Altitude  =  19.7  km 
(64.5  kft) 

Altitude  =  35.1  km 
(115  kft) 

Flight  Mach  Number 

5.0 

5.0 

Static  Pressure,  Pa 
(atm) 

5833.7  (0.57) 

558.2  Pa  (0.055) 

Air  Density,  kg/m° 

0.0938  kg/m3  ' 

0.0082  kg/mJ 

Static  Temperature,  K 

216.67 

237.2 

Missile  Velocity,  m/sec1 

1475.0 

1543.3 

Angle  of  Attack,  deg 

0.0 

0.0 

Yaw  Angle,  deg 

0.0 

0.0 

Molecular  Weight 

28.967 

28.967 

Ratio  of  Specific  Heats 

1.4 

1.4 

the  center  of  gravity  on  the  interceptor  missile  and 
discharging  perpendicular  to  the  Mach  5  free- 
stream  flow.  The  CFD  investigation  included  four 
solutions  at  each  of  the  altitudes  mentioned  above. 
The  solutions  comprised  steady-state  calculations 
with  the  divert  thruster  turned  jet  off  and  a  steady- 
state  solution  with  the  jet  on,  as  well  as  the  tran¬ 
sient  startup  of  the  jet  and  a  transient  shutdown  of 
the  jet.  The  eight  CFD  cases  investigated  are  sum¬ 
marized  in  Table  2. 


Table  2.  Jet  Interaction  CFD  Cases 


Case 

No. 

Time 

Integration 

Jet  Flow 

Mach 

No. 

Altitude  kft 
(km) 

1 

Steady  State 

Off 

5 

64.5  (19.7) 

2 

Steady  State 

On 

5 

64.5  (19.7) 

3 

Transient 

Jet  Turning  On 

5 

64.5  (19.7) 

4 

Transient 

Jet  Turning  Off 

5 

64.5  (19.7) 

5 

Steady  State 

Off 

5 

115.0(35.1) 

6 

Steady  State 

On 

5 

115.0  (35.1) 

mi 

Transient 

Jet  Turning  On 

5 

115.0  (35.1) 

HZ 

Transient 

Jet  Turning  Off 

5 

115.0(35.1) 

Geometry  and  Computational  Grid 


The  generic  interceptor  missile  geometry  was 
specified  as  a  relatively  long  and  slender  configura¬ 
tion  (L/D  =  14.1)  equipped  with  fixed  dorsal  and  tail 
fins.  The  divert  thruster  nozzle  exit  was  positioned 
flush  with  the  missile  surface  between  the  dorsal 
fins  at  the  center  of  gravity  of  the  missile.  The  total 
length  of  the  interceptor  missile  was  specified  as 
195.74  in.,  which  is  equivalent  to  14.5  calibers.  One 
caliber  is  equivalent  to  one  missile  diameter,  speci¬ 
fied  as  13.5  in.  (34.3  cm).  The  nose  geometry  was 
modeled  as  a  2  caliber  tangent  ogive  with  a  0.33-in. 
(0.84  cm)  radius  spherical  nose  cap.  The  afterbody 
section  is  a  12.5-caliber  cylinder.  Very  thin  dorsal 
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and  tail  fins  are  rectangular  in  cross  section  and 
symmetrically  positioned  around  the  body  at  45  and 
135  deg  from  the  symmetry  plane.  The  dorsal  fins 
have  a  5-caliber  chord  and  a  0.5-caliber  span  with 
the  leading  edge  located  5  calibers  from  the  nose 
tip.  The  tail  fins  have  a  1 -caliber  chord  and  a  0.5- 
caliber  span,  with  the  leading  edge  located  13.5 
calibers  form  the  nose  tip.  The  circular  divert  jet 
thrust  chamber  exit  slot  is  flush  to  the  missile  body 
at  90  degrees.  The  exit  slot  diameter  is  2.964  in. 
(7.53  cm).  As  shown  in  Fig.  1 ,  the  center  of  the  exit 
slot  is  located  on  the  interceptor  missile  body  6  cal¬ 
ibers  form  the  missile  nose  tip  in  the  vehicle  sym¬ 
metry  plane  between  the  dorsal  fins. 

A  computational  grid  was  constructed  for  this 
geometry,  including  the  interceptor  missile  and  the 
internal  thrust  chamber  of  the  divert  thruster.  The 
grid  consists  of  13  blocked  meshes  (approximately 
2.6  million  grid  points)  for  the  low-altitude  cases 
and  14  blocked  meshes  (approximately  3  million 
grid  points)  for  the  high-altitude  cases.  An  addi¬ 
tional  mesh  was  required  for  the  high-altitude 
cases  to  completely  capture  the  relatively  larger 
underexpanded  jet  plume  that  occurs  at  the  higher 
altitude,  lower  ambient  pressure  environment.  The 
computational  surface  grid  for  the  interceptor  mis¬ 
sile  is  presented  in  Fig.  2.  A  lateral  control  jet,  also 
known  as  a  divert  control  jet,  is  located  at  the  mis¬ 
sile  center  of  gravity  between  the  dorsal  fins.  The 
jet  thrust  is  produced  by  a  convergent-divergent 


Fig.  2.  Surface  grid. 


nozzle  having  supersonic  exit  conditions.  The  nom¬ 
inal  exit  Mach  number  is  3.0.  An  enlargement  of 
the  grid  near  the  divert  jet  slot  location  on  the  mis¬ 
sile  surface  is  shown  in  Fig.  3.  A  very  high  grid 
density  was  constructed  in  the  region  upstream  of 
the  divert  jet  exit  to  obtain  spatial  details  of  the  flow 
gradients  in  this  complex  Jl  region.  A  pitch  plane- 
of-symmetry  was  assumed  for  all  cases.  The  vol¬ 
ume  grid  for  the  free-stream  computational  domain 
including  all  blocked  meshes,  is  shown  in  Fig.  4. 
The  divert  jet  thruster  internal  flow  field  was  com¬ 
puted  using  the  grid  shown  in  Fig.  5.  Grid  points 
were  packed  near  the  nozzle  walls  to  capture  vis¬ 
cous  effects.  The  thruster  flow  field  was  simulated 


from  specified  total  conditions  in  the  chamber, 
assuming  thermally  perfect,  y  =  1.4  flow.  This 
enabled  the  computations  to  capture  the  transient 
effects  of  the  jet  startup  and  account  for  the  jet 
boundary  layer  and  divergence  at  the  exit  plane. 


Fig.  3.  Closeup  view  of  surface  grid  around  nozzle 
exit. 


Fig.  4.  Volume  grid  with  high-altitude  auxiliary  grid. 


Fig.  5.  Volume  grid  for  lateral  thruster. 
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Analysis 

The  General  Propulsion  Analysis  Chemical 
Kinetic  and  Two-Phase,  GPACT3  was  applied 
exclusively  to  simulate  the  complex  at  each  condi¬ 
tion.  GPACT  solves  the  integral  form  of  the  three- 
dimensional,  time-dependent,  Reynolds-averaged, 
full,  compressible  Navier-Stokes  equations.  GPACT 
is  capable  of  solving  subsets  of  these  equations, 
including  two-dimensional  and  axisymmetric  equa¬ 
tions.  GPACT  includes  an  optional  multizone  capa¬ 
bility,  which  was  applied  in  this  study.  The  algorithm 
employs  a  shock-capturing  solution  scheme.  This 
methodology  is  fully  conservative  because  of  the 
consistent  manner  in  which  properties  such  as  com¬ 
putational  cell  volumes,  surface  areas,  and  numeri¬ 
cal-flux  functions  are  evaluated.  The  flow  solver 
contains  thermodynamic  models  and  chemistry 
models.  High-temperature  effects  such  as  vibra¬ 
tional  relaxation  and  preferential  dissociation  can  be 
included  in  the  calculations.  For  this  investigation,  a 
thermally  perfect  gas  was  assumed.  Several  turbu¬ 
lence  models  are  provided  as  user-selected  options. 
The  turbulence  models  include  the  two-layer  alge¬ 
braic  model  of  Baldwin  and  Lomax,  and  fully  cou¬ 
pled,  two-equation  k-e  models.  An  approximation 
accounting  for  turbulence-chemistry  interacts  is  also 
available.  The  current  study  uses  Van-Leer  flux 
splitting,  assumes  the  flow  medium  is  a  thermally 
ideal  gas  (y  =  1 .4),  and  a  k-e  model.  Forces  and 
moments  on  the  missile  body  were  obtained  by  inte¬ 
grating  the  surface  pressures  and  viscous  shear 
stresses  over  the  entire  surface  of  the  missile.  The 
effects  of  the  interceptor  base  entrainment  and 
associated  forces  were  not  evaluated.  Divert  jet 
thrust  was  obtained  by  integrating  the  calculated 
flux  across  the  nozzle  exit  diameter  of  the  jet. 

Results  and  Discussion 

The  GPACT  model  provided  converged  flow- 
field  solutions  for  all  cases  considered  in  the  numer¬ 
ical  investigation.  Computational  times  for  con¬ 
verged  steady-state  solutions  ranged  from  100-150 
CPU  hours  on  a  Hewlett  Packard  SPP-2000  Exem¬ 
plar  computer  using  eight  processors.  Less  CPU 
time  was  required  for  solutions  obtained  when  the 
divert  thruster  was  not  operating.  The  transient 
computations  were  executed  in  time-accurate  mode 
and  required  considerably  longer  CPU  times.  Addi¬ 
tionally,  transient  CFD  analyses  cannot  take  advan¬ 


tage  of  convergence  acceleration  techniques.  The 
transient  computations  required  20-40  CPU  hours 
to  integrate  over  1  msec  of  physical  time,  depend¬ 
ing  on  the  time  step  and  grid  resolution.  Initially,  a 
time  step  of  0.2  psec  was  used  for  the  transient 
solutions.  This  value  was  increased  in  a  stepwise 
fashion  to  0.5  psec  as  the  computation  proceeded. 
To  assure  an  eventual  steady-state  condition,  the 
transient  methodology  was  continued  to  14  psec. 

Free-stream  conditions  were  specified  on  the 
far-field  boundaries,  and  viscous  no-slip  assump¬ 
tions  were  specified  on  all  missile  surfaces,  includ¬ 
ing  the  divert  jet  thrust  chamber  walls.  The  thrust 
chamber  wall  boundary  conditions  specified  for  the 
divert  nozzle  inflow  depended  on  the  whether  the 
thruster  jet  was  operational.  In  the  jet-off  case,  the 
nozzle  inflow  boundary  was  replaced  with  a  slip 
wall,  otherwise;  a  subsonic  inflow  boundary  condi¬ 
tion  was  used  for  the  jet-on  simulations.  The  divert 
thruster  stagnation  conditions  (P0j  =  2590  psia  and 
Toj  =  4100°R)  were  specified  for  the  jet  subsonic 
inflow  boundary  condition.  Outflow  boundary  condi¬ 
tions  equal  to  the  free-stream  ambient  conditions 
were  specified  for  the  downstream  boundary.  For 
those  simulations  where  the  divert  jet  was  turned 
off,  the  empty  thrust  chamber  cavity,  embedded 
within  the  interceptor  surface,  was  included  in  the 
computational  domain.  The  interceptor  missile  base 
aerodynamic  effects  were  not  computed;  therefore, 
the  force  and  moment  contributions  of  the  base 
region  are  neglected. 

Flow-field  results  from  the  steady-state,  jet-on 
calculation  are  shown  in  Fig.  6.  Mach  number  con¬ 
tours  at  radial  cross  section  are  displayed  at  three 
axial  locations  flow  corresponding  to  axial  positions 
on  the  missile  body:  (1)  just  downstream  of  the 
divert  jet,  (2)  the  end  of  the  dorsal  fins,  and  (3)  the 
end  of  the  missile.  The  behavior  of  the  under¬ 
expanded  divert  jet  is  indicated  by  the  rapid  diver¬ 
gence  of  the  jet  exhaust  in  the  first  cross  section. 
Effects  of  the  divert  jet  and  free-stream  interactions 
are  noticeable  at  the  downstream  locations.  Note 
the  high-pressure  interaction  region  impinging  on 
the  fins.  This  integrity  of  flow-field  results  in  high  gra¬ 
dient  regions,  such  as  the  Jl  zones,  is  strongly 
dependent  on  the  computational  grid.  Including  the 
divert  thrust  chamber  in  the  calculation,  all  the 
effects  of  the  jet  boundary  layer  are  simulated. 
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Fig.  6.  Mach  number  contours  for  64.5-kft  altitude, 
jet-on  case. 
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Since  the  divert  jet  operates  in  an  underex¬ 
panded  mode  at  both  altitude  conditions,  the  flow  at 
the  nozzle  exit  plane  is  supersonic  (approximately 
Mach  number  3.0).  Under  these  conditions,  the 
internal  thrust  chamber  flow-field  conditions  are 
independent  of  the  free  stream,  resulting  in  identi¬ 
cal  steady-state  internal  thrust  chamber  flow  fields 
at  both  altitudes.  However,  the  thrust  at  the  higher 
altitude  is  expected  to  be  larger  because  of  the 
smaller  ambient  force  term  at  the  nozzle  exit.  Com¬ 
puted  velocity  vectors  from  the  jet-on  steady-state 
solution,  including  the  thrust  chamber  and  immedi¬ 
ate  free-stream  interaction  regions  for  the  lower  alti¬ 
tude  (64.5  kft)  Mach  number  5  free-stream  condi¬ 
tions  are  shown  in  Fig.  7.  The  large  circulation  bub¬ 
ble  upstream  created  by  the  free-stream  interaction 
with  the  divert  jet  plume  is  evident  in  the  region 
upstream  of  the  divert  nozzle  exit  slot.  A 
secondary  separation  region  is  also 
apparent.  These  recirculation  zones  cre¬ 
ate  a  high-pressure  region  and  an  associ¬ 
ated  force  on  the  missile  surface.  The 
force  exerted  by  the  recirculation  effect 
magnifies  the  divert  jet’s  aerodynamic 
force,  and  also  influences  the  moment 
about  the  center  of  gravity  of  the  intercep¬ 
tor  missile.  As  indicated  on  Fig.  7,  the 
interaction  on  the  downstream  side  of  the 
divert  jet  creates  an  expansion  region.  Fig¬ 
ure  8  shows  a  comparison  of  64.5-  and 
115-kft  Mach  number  contours.  The  lower 
ambient  pressure  at  the  high  altitude  pro¬ 
duces  a  significantly  larger  jet  plume  and 
separation  region  compared  to  the  lower- 
altitude  results. 


Fig.  7.  Velocity  vectors  for  64.5-kft  altitude,  jet  on 
case. 


Fig.  8.  Jet  interation  Mach  No.  contours-comparison 
of  64.5- and  115-kft. 


Fig.  9.  Surface  pressure  at  64.5-ktf  altitude,  jet  turning  on  case. 
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Fig.  10.  Mach  number  contours  for  64.5-ktf  attitude, 
on  case. 

Figures  9-16  are  the  results  of  the  transient  Jl 
solution  at  the  64.5-kft  altitude  condition  showing 
the  effects  of  the  divert  jet  startup  process.  Several 
transient  time  slices  are  shown  culminating  in 
steady-state  flow  conditions  at  approximately  6 
psec.  Figure  9  depicts  the  transient  time  slices  of 
calculated  static  pressure  on  the  surface  of  the 
interceptor  vehicle  from  initiation  of  the  jet  flow  at 
time  =  0.0  to  established  steady-state  flow  at  8 
psec.  Note  the  high-pressure  interaction  region 
impinging  on  the  fins  at  3.0  seconds.  Conse¬ 
quences  of  this  impingement  will  be  observed  in 
subsequent  force  analysis.  Figure  10  includes  cor¬ 
responding  Mach  number  contours  in  the  free- 
stream  and  Jl  region  symmetry  plane  covering  the 
same  time  span  as  the  flow  develops.  Steady-state 
Mach  number  contours  are  displayed  at  8.0  psec. 
These  results  clearly  shown  the  development  of  Jl 
shock  and  the  associated  high  surface  pressures 
and  low  Mach  numbers  in  the  recirculation  zone 
upstream  of  the  divert  jet  slot.  The  downstream 
expansion  regions  are  characterized  by  supersonic 
flow  and  decreasing  pressure.  It  is  clear  that  Mach 
number  contours  on  the  symmetry  plane  are  shown 
in  Fig.  10  over  the  period  of  the  jet  turning  on.  A 
small  disturbance  can  be  seen  at  t  =  0,  caused  by 
the  presence  of  the  nozzle  cavity  with  no  flow. 

The  surface  integrated  normal  force  and  pitch¬ 
ing  moments  calculated  for  the  transient  jet  startup 


|  process  (64.5  kft)  are  shown  in  Figs.  1 1 
and  12,  respectively.  At  t  =  0,  the  jet  is 
]  turned  on  and  the  force  rises  abruptly 
due  to  the  jet  force  alone  (-920  N).  The 
exerted  force  continues  to  rise  as  the 
j  high-pressure  recirculation  region  devel¬ 
ops  and  spans  increasing  areas  on  the 
cylindrical  missile  body.  The  enhanced 
transient  effect  begins  at  approximately 
2.6  seconds  and  occurs  over  a  1 .2-psec 
burst,  creating  normal  body  forces  in 
excess  of  3300  N  at  the  maximum  point. 
A  steady-state  condition  is  established 
near  6  psec  The  steady-state  force  is 
approximately  1750  N  and  is  greater 
than  the  magnitude  of  the  jet  force  alone 
(920  N).  The  pitching  moment  (computed 
jet  turning  about  the  center  of  the  divert  jet)  quickly 
rises  to  a  value  of  550  N-m  as  the  jet  is 
turned  on,  and  then  appears  to  remain 
steady  over  the  period  from  0-2.8  msec.  The  pitch¬ 
ing  moment  then  rises  to  a  value  approaching  2800 
N-m  over  the  4-sec  interval  from  3-6.8  psec.  At 
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Fig.  1 1 .  Integrated  normal  force  for  jet  turning  on  at 
64.5  kft. 
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Fig.  12.  Integrated  pitching  moment  for  jet  turning 
on  at  64.5  kft. 
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approximately  3.8  msec,  the  jet  interaction  region 
impinges  on  the  dorsal  fins  in  an  abrupt  manner 
with  a  force  centered  at  the  axial  location  of  the  jet. 
This  is  evident,  as  there  is  an  abrupt  change  in 
moment  with  little  change  in  force.  After  4.0  msec, 
the  high-pressure  region  established  in  front  of  the 
divert  jet  continues  to  travel  upstream,  causing  a 
slight  drop  in  surface  pressure  as  the  moving 
shocks  begin  to  settle  down  and  become  more  sta¬ 
tionary  and  oblique  to  the  oncoming  free-stream 
flow.  This  causes  a  net  decrease  in  the  magnifica¬ 
tion  force  caused  by  the  separation  region  acting 
on  the  missile.  However,  the  resulting  steady-state 
force  (1750  N)  is  greater  than  double  the  force 
value  computed  for  the  divert  jet  alone  (920  N). 

Figures  1 3  and  1 4  show  the  integrated  force  and 
moment  for  the  jet  shutdown  process  at  the  64.5-kft 
conditions.  Time  equal  to  0.0  corresponds  to  when 
the  thruster  jet  flow  is  extinguished.  The  force 
immediately  begins  to  decrease  as  the  jet  shuts 
down.  The  interaction  force  appears  to  remain 
steady  for  1.1  msec  as  the  flow  field  adjusts  to  the 


new  condition.  The  moment  also  appears  to  experi¬ 
ence  a  delayed  response  to  the  jet  shutdown,  and 
begins  to  decrease  after  a  period  of  approximately 
1.5  msec.  The  final  steady-state  force  equal  to 
approximately  80  N  and  moment,  equal  to  120  N-m, 
are  the  result  of  the  steady  free  steam  flow.  The 
normal  force  and  pitching  moment  over  a  complete 
cycle,  from  jet  start  to  shutdown  at  64.5  kft,  are 
shown  in  Figs.  15-16,  respectively. 
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Fig.  15.  Integrated  normal  forces  from  start  to  shut¬ 
down  at  64.5  kft. 
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Fig.  13.  Integrated  normal  force  for  jet  turning  off  at 
64.5  kft. 
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Fig.  16.  Integrated  pitching  moment  time  trace  from 
start  to  steady-state  shutdown. 
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Fig.  14.  Integrated  pitching  moment  for  jet  turning 
off  at  64.5  kft. 


Corresponding  results  for  the  transient  high-alti¬ 
tude  solutions  are  shown  in  Figs.  17-21.  The 
results  have  very  similar  trends  when  contrasted 
with  the  low-altitude  results.  Steady-state  velocity 
vectors  for  the  jet-on  calculation,  showing  the 
divert  jet  thrust  chamber  and  immediate  free- 
stream  regime,  is  shown  in  Fig.  17.  Since  the  same 
jet  conditions  were  used  for  this  higher-altitude 
case,  the  divert  jet  nozzle  is  even  more  underex¬ 
panded  at  this  higher  altitude  condition,  resulting  in 
a  larger  exhaust  plume  and  deeper  jet  penetration 
into  the  free  stream.  The  higher  altitude  case 
exhibits  a  larger  separation  region  because  of  the 
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reduced  momentum  of  the  free  stream,  providing 
conditions  that  are  more  favorable  for  separation. 
Comparison  of  higher  altitude  at  1 1 5  kft  with  lower 
altitude  at  64.5  kft  are  shown  in  Fig.  8.  Note  that 
although  the  separation  region  is  increased  at  this 
altitude,  the  resulting  pressure  increase  on  the  mis¬ 
sile  body  is  lower  when  compared  with  the  lower 
altitude  case.  As  will  be  shown  in  subsequent  force 
profiles,  the  lower  pressures  in  the  separation  zone 
(compared  with  the  64.5-kft  solution)  results  in  a 
decreased  force  magnification  effect. 
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Figures  18  and  1 9  show  the  transient  jet  startup 
process  for  the  high-altitude  case.  Static 


Fig.  17.  Velocity  vectors  for  115-kft  altitude,  jet  on 
case. 


missile  surface  pressures  are  shown  in 
Fig.  18  at  time  slices  corresponding  to 
initial  startup  until  steady-state  flow  is 
reached  at  approximately  8.0  psec.  The 
plume  is  bigger  at  this  higher  altitude 
condition,  and  the  Jl  region  impinges  on 
the  dorsal  fins  earlier.  However, 
because  of  the  lower  pressure  environ¬ 
ment,  it  is  expected  that  the  correspond¬ 
ing  impingement  force  will  be  lower. 
Free-stream  Mach  number  contours  in 
the  horizontal  symmetry  plane  of  the 
divert  nozzle  over  the  same  time  span 
are  shown  in  Fig.  19.  Larger  plume 
expansion,  deep  divert  jet  penetration 
into  the  free  stream,  and  very  large  recir¬ 
culation  regions  are  evident  in  these 
results. 


The  surface  integrated  force  and 
moment  values  for  the  jet  startup  tran¬ 
sient  process  at  the  higher  altitude  (115 
kft)  conditions  are  shown  in  Figs.  20-21 . 
These  plots  exhibit  similar  trends  to  the 
lower  altitude  results,  but  are  less  abrupt. 
The  transient  spike  occurs  earlier  in  the 
startup  process  (3.5  msec  vs.  3.8  msec) 
and  is  more  gradual  compared  with  the 
lower-altitude  solution.  The  Jl  impinge¬ 
ment  on  the  dorsal  fins  occurred  earlier  at 
this  condition;  the  transient  burst  is  pre¬ 
dicted  to  occur  earlier.  Recall  that  the 
time  corresponding  to  the  peak  force  was 
correlated  with  the  impingement  of  the  Jl 
region  on  the  missile  dorsal  fines  at  the 
lower  altitude  condition. 


Fig.  19.  Mach  number  contours  for  115-kft  altitude,  jet  turning  on 
case. 
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Fig.  20.  Integrated  normal  force  for  jet  turning  on  at 
115  kft. 
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Fig.  21.  Integrated  pitching  moment  for  jet  turning 
on  at  115-kft. 

The  peak  transient  interaction  force  (-2200  N) 
at  the  higher-altitude  conditions  is  considerably 
less  than  the  lower-altitude  peak  force  (3300  N). 
The  is  attributed  to  the  lower  pressures  at  the 
higher  altitude.  The  moment  transient  increase  is 
very  gradual,  lagging  the  force  trend,  as  was  evi¬ 
dent  in  the  lower-altitude  results. 

Force  and  moment  profiles  for  the  jet  shutdown 
transient  at  1 15  kft  altitude  is  shown  in  Figs.  22  and 
23,  respectively.  These  results  are  very  similar  to 
the  shutdown  behavior  at  the  lower  altitude.  The 
steady-state  value  at  the  end  of  the  shutdown  tran¬ 
sient  is  lower  because  of  the  lower  pressure  in  the 
free  stream.  These  results  indicate  that  shorter 
times  are  required  to  reach  steady-state  conditions 
at  the  high-altitude  condition.  This  could  be  influ¬ 
enced  by  the  very  large  jet  penetration  and  recircu¬ 
lation  regions  that  result  at  higher  altitudes.  The 
normal  force  andpitching  moment  over  a  complete 
cycle  from  jet  start  to  shutdown  at  115  kft  are 
shown  in  Figs.  24  and  25,  respectively. 
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Fig.  22.  Integrated  normal  force  for  jet  turning  off  at 
115-kft. 
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Fig.  23.  Integrated  pitching  moment  for  jet  turning 
off  at  115  kft. 
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Fig.  24.  Integrated  normal  forces  from  start  to  shut¬ 
down  at  115  kft. 
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Fig.  25.  Integrated  pitching  moment  from  start  to 
shutdown  at  115  kft. 


10 

UNCLASSIFIED 


UNCLASSIFIED 


Summary 

Transient  jet  interaction  computations  were 
completed  to  assess  the  transient  effects  of  a  reac¬ 
tion  control  thruster  on  the  free-stream  flow 
field, and  normal  forces  and  pitching  moments  on  a 
generic  interceptor  missile.  The  results  indicate  a 
strong  transient  effect  that  must  be  considered 
when  designing  control  algorithms  for  pulsed-jet 
reaction  control  systems  with  short  pulse  times. 
The  transient  interaction  at  low  altitudes  on  this 
type  of  finned  geometry  causes  short  durations  of 
high  force  that  diminish  at  higher  altitudes.  This 
very  strong  jolt  is  caused  by  the  high-pressure 
region  in  front  of  the  jet  exhaust  gasses  as  it 
abruptly  impacts  the  dorsal  fins.  These  pressures 
are  higher  than  the  steady  pressures  because  of 
moving  shocks.  The  transient  effects  seem  to  last 
approximately  3-6  msec  at  both  high  and  low  alti¬ 
tudes  before  the  force  begins  to  steady  out.  The  jet 
shutdown  case  required  less  time  to  reach  steady- 
state  conditions  than  the  jet  startup  process. 

This  investigation  demonstrates  the  utility  of 
state-of-the-art  CFD  tools  to  address  issues  per¬ 
taining  to  the  most  complex  flow  interactions  of 
interest  in  the  propulsion  arena.  While  these  solu¬ 
tions  require  considerable  CPU  resources,  they 
provide  an  economical  means  to  supplement  inter¬ 
ceptor  design  and  testing  and  have  the  potential  to 
save  development  and  testing  resources.  CFD 
solutions  such  as  those  presented  in  this  study  can 
be  applied  to  maximize  the  test  matrices  and  to 
evaluate  paper  designs  and  identify  potential  prob¬ 
lem  area  prior  to  production.  The  CFD  flow  fields 
can  be  analyzed  to  gain  insight  into  the  force 
amplication  phenomena  and  to  understand  the 
transient  behavior.  The  flow-field  details  calculated 
and  visualized  in  this  study  cannot  be  measured  at 
the  harsh  conditions  using  modern  flow  visualiza¬ 
tion  techniques.  Additional  CFD  studies  investigat¬ 
ing  the  effects  of  chemistry  in  the  interaction  region 
are  needed  to  complete  the  CFD  feasibility  demon¬ 
stration. 
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